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Abstract. Using an approach that treats the Ricci scalar 
CN ■ itself as a degree of freedom, we analyze the cosmological 

evolution within an f(R) model that has been proposed 
recently (exponential gravity) and that can be viable for 
explaining the accelerated expansion and other features of the 
1 Universe. This approach differs from the usual scalar-tensor 

. method and, among other things, it spares us from dealing 

with unnecessary discussions about frames. It also leads to a 
simple system of equations which is particularly suited for a 
numerical analysis. 



1. Introduction 

Since the discovery of the accelerated expansion of the Universe from the 
observations of supernovae la [T7J [TBI d] an d its interpretation using the 
ACDM model of standard cosmology, a large amount of investigations 
have been devoted to explain the same phenomenon but using dark energy 
substances different from the cosmological constant A. One can rank the 
dark energy models from the most "conservative" to the most "radical" 
ones. Among the former we can mention those which do not introduce 
new fields or modifications to general relativity but which consider that 
inhomogeneities in the Universe could be enough to account for the 
observations [TS]. There are also the models that introduce new fields 



and perfect fluids with exotic equations of state within the framework of 
general relativity (GR) in order to avoid the "problems" associated with 
A |19[ 114) For instance, quintessence, k-essence, and Chaplygin 
gas are some of the most popular models of this kind. The most radical 
attempts to explain the accelerated expansion are perhaps those which 
propose to modify GR while keeping the hypothesis of homogeneity and 
isotropy as a first approximation. Many alternate theories of gravity 
have been proposed to explain this phenomenon as well as those related 
with the dark matter (e.g. rotation curves of galaxies). The modified 
f(R) metric gravity is just one of such theories and maybe the most 
analyzed one in the last ten years, where the geometry takes care of 
mimicking the dark energy. This alternative is certainly radical since GR 
has been thoroughly tested for almost one hundred years and it has not 
only supported all the tests but in addition most of its predictions have 
been confirmed as well. Thus, the challenge of modified gravity is both to 
be consistent with GR tests and also to explain the phenomena they were 
called for. This is a no trivial task and many f(R) models have failed in 
the attempt. The story concerning dark substances does not end with the 
accelerated expansion. The measurements of the angular distribution of 
cosmic background radiation anisotropies in the sky can also be explained 
by the ACDM model, and therefore, the task for the alternative models, 
theories or dark energy substances is even more demanding. 

As we mentioned, f(R) theories have been studied in detail in a recent 
past and it is out of the scope of the present article to discuss all the 
properties, problems and features associated with some of the specific 
models proposed before (see Refs. [T6j [101 EOl SI [6] for a review). 

Our aim is to report the results of a potentially viable candidate, termed 
exponential gravity, as a model for the accelerated expansion, but using an 
approach that has been proposed recently by us [TT] and which avoids 
the identification with the scalar-tensor theories. The reason to follow 
this "unorthodox" method is because in some cases the scalar-tensor (ST) 
method can lead to ill-defined potentials, and moreover because we want to 
circumvent any possible discussion concerning the use of frames (Einstein 
vs Jordan). Debates of this sort plague the subject, some of which have 
only led to create confusion instead of shedding light. 

With our technique we propose to treat the Ricci scalar itself as a 
degree of freedom, instead of using (j> — Jr as in the ST method (hereafter 
a subindex R indicates d/dp). Our approach also spare us of inverting 
all the quantities depending on R for treating them as functions of </>. 
Moreover, we have found that in several specific applications the field 
equations can be recasted in a rather friendly way that allows us to 
treat them numerically or even analytically [12] E] • In the next section 
we present our method and apply it to the Friedmann-Roberson- Walker 



(FRW) spacetime within the scope of analyzing the cosmological evolution 
using the exponential gravity model. The analysis of other viable f{R) 
models using the current approach can be seen in and in references 
therein using other techniques. 

2. f(R) theories, the Ricci scalar approach 

The action in f(R) gravity is given by: 

S[g ab ^}^J^V^9d 4 x + S mM [gab,'4'} , (1) 

where n = 8ttGq (we use units where c = 1), f(R) is a sufficiently 
diffcrcntiable but otherwise a priori arbitrary function of the Ricci scalar 
R, and %p represents schematically the matter fields. The field equation 
obtained from Eq. (fTJ) is: 

fnRab ~ T^fgab - (V a Vb ~ 9ab^) f R, = nT ab , (2) 

where f R indicates d R f, □ = g ab V a V b is the covariant D'Alambertian and 
T a b is the energy-momentum tensor (EMT) of matter associated with the 
ip fields. From this equation it is straightforward to obtain the following 
equation and its trace |12[ [TT] 



G a b — ~f 

JR. 



fRRV a V b R+f R R R (y a R)(V b R)-^[],'h,+j 2,,T')h-/,-T„,., 



DR : 



(3) 

K T-3f RRR (VR) 2 + 2f-Rf R , (4) 



SfRR 

where (Vi?) 2 := g ab {V a R){V b R) and T := T a a . 

The idea is then to solve simultaneously Eqs. ^ and d3| for the metric 
g ab and R as a system of coupled partial differential equations. 

It is important to mention that the field equations imply that the EMT 
of matter alone is conserved, i.e., it satisfies V Q T ab = 0. 

In this contribution we shall focus on the model f(R) = — 
A(l — e~^)], referred to as exponential gravity, where R = A is 

a positive dimensionless constant and i?* > is also a constant that 
fixes the built-in scale and which is of the order of the current Hubble 
parameter Hq. This kind of exponential models have been analyzed in 
the past by several authors using a different technique [5] Q~31 HU [H [HI H] . 
Other variants of this model have also been analyzed [21]. The scalar 
f R = l — \e~ R is positive provided R > i?*lnA and this condition ensures 
G ff := Go/f R > 0. This latter is always satisfied in the cosmological 
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Figure 1. (color online) f(R) exponential gravity (left panel) and the 
potential V(R) (right panel) for several values of A. 



solutions given below. The possible de Sitter points correspond to trivial 
solutions R = Ri = const, of Eq. (j4| in vacuum (T = 0) and give rise to an 
effective cosmological constant A c ff = i?i/4 in Eq. ^ {G a b = —gables in 
vacuum). Here R\ > is a critical point of the "potential" V(R) such that 
Vr{R!) = with Vr := (2f-RJ R )/3 = R*[R(l + \e- R )-2\(l-e- k )}/3, 
and Vrr = 1 - A(l + R)e- R . The "potential" is given by V(R) = 
Rl[R{R - 4A) - 2\(R + 3)e~ R )}/6. For < A < 1 one can easily see that 
V(R) has just one critical point at R — which is not a de Sitter point as 
in this case A e g = 0. The point is a global minimum (c.f. Vrr(0) = 1 — A). 
For A > 1 there is a local maximum at R = and a global minimum 
at R = Ri > which corresponds to the actual de Sitter point that the 
cosmological solution reaches asymptotically in the future. There is also 
a local minimum at R = i?2 < 0, but it is an anti de Sitter point which is 
never reached as the cosmological solutions take place only in the domain 
R > 0. The potential V(R) is depicted in Figure [1] (right panel) where 
one can appreciate the critical points just described. In the high curvature 
regime R ^> 1, we have f(R) ~ R — Ai?*, and thus the model acquires an 
effective cosmological constant A£g := Ai?*/2 (c.f. left panel of Fig. [TJ. 
From the figure we see that for A sufficiently high, the de Sitter point 
verifies Ri 3> 1, and thus R\ w 2Ai?„ as it turns out from Vr(Ri) = 0. 
Therefore A e g rj A^. Finally, we stress that Jrr — \e~ R /R* > which 
ensures that no singularities are found in the equations due to this scalar 
and moreover it guarantees that given Vrr(Ri) > the effective mass 
to 2 := Vrr/Srr = (,f R - RfRR)/(3f RR )\ Rl of the scalar mode is positive. 



3. Cosmology in f(R) 



We assume the spatially flat FRW metric given by: 

ds 2 = -dt 2 + a 2 {t) [dr 2 + r 2 (dO 2 + sin 2 ddip 2 )] . (5) 
From Eqs. §S§ and (j3J) we have, 



R = -3HR - ' 



3f RRR R 2 + 2f-f R R + KT , (6) 



H 2 = ^( P + PX ) , (7) 



if = ~H 2 - 



£ {p + PX + 3(prad + Px)} ■ (8) 



where a dot stands for d/dt and H = a/ a, is the Hubble expansion. In 
the above equations we have included the energy density p associated with 
matter (baryons, dark matter and radiation) as well as the GDE density 
px and pressure px given respectively by [IT] 

px = { \ {IrR ~ f) ~ 3fRiiHk +k p^- } > ( 9 ) 

Px = - 7V ( k (f* R + f) + ZfunHR - k (p - 3p rad / R ) ) . (10) 

Another differential equation that can be used to solve for if instead 
of Eq. © is given by R = 6 (if + 2if 2 ). This latter is no other than 
the Ricci scalar computed directly from the metric (j^J. Equation Q 
amounts to the modified Hamiltonian constraint which we use to set the 
initial data and also to monitor the accuracy of the numerical solutions 
at every integration step. At this regard, we stress that we shall not use 
the cosmic time t but instead a = ln(a/ao) as "time" parameter (see 
Ref. [H]), where oo is the present value of a. Notice that at the de Sitter 
point R — > Ri = const, where 2/(i?i) = Rif R (Ri) and with p — > 0, 
Prad — » Eqs. © and (fTU)) . lead to px — > A e ff/« and px — > —A e g/K, 
respectively, and from Eqs. © and ©, if 2 -> if 2 ac = A cff /3 = fli/12, 
and g = —a/{aH 2 ) — > q vac = —1. So the main idea behind all /(f?) models 
is that as the Universe evolves, R — > i?i, and thus the GDE dominates 
and mimics an effective cosmological constant that allows to explain the 
accelerated expansion required to account for the observations. 

The matter variables obey the conservation equation pi = 
—3H(pi + Pi) for each fluid component (with pbar,DM = and p ra( j = 
Prad/3) which integrates straightforwardly and gives rise to the usual ex- 
pression for the energy density of matter plus radiation: p — (Pb ar + 
/9p M )(a/ao)~ 3 + Prad( a / a o)~ 4 i where the knotted quantities indicate their 



Figure 2. (color online) Ricci scalar (left panel) and the Hubble expansion 
(right panel) for several values of A and i?* (given in units of Hq). 
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Figure 3. (color online) Evolution of fi ma tt (red solid line) and tlx (blue 
solid line) for A = 1 and i?* = AH§ (left panel) and A = 2.5 and R* = 2Hq 
(right panel). For reference the corresponding quantities of the ACDM 
model are included in each panel (dashed lines). 



values today. The X-fluid variables ([H]) and ([TU]) also satisfy a conserva- 
tion equation similar to the one above, but with an EOS uix '■= Px/px 
that evolves in cosmic time. Other possible inequivalent definitions of 
Px, Px and lox have been adopted in the past, but they suffer of several 
drawbacks (see [TT] for a detailed discussion). 

The total EOS is defined by u> to t = (ftad + Px)/(p + Px) which using 
Eqs. dU) and {ID} yields 

IUrR- f)-3f RR HR + Kp\ 

This EOS allows us to track the epochs where the Universe is expanding 
in a decelerating or accelerating fashion. If io to t < — 1/3 then a > 0, while 
a < if Wtot > —1/3. 
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Figure 4. (color online) The EOS uix (left panel) and the total EOS Wtot 
(right panel). 

4. Numerical Results and Discussion 

We integrate the differential equations forward, from past to future, 
starting from a given z = do/a — 1, where matter dominates, to z — > — 1 
where the GDE prevails. The initial conditions are fixed as in [IT]. Figure^] 
shows the Hubble expansion and the Ricci scalar for several values of A. 
In all the cases where A > 1, R reaches the de Sitter point at the global 
minimum of V(R). For A = 1 because the potential is very flat around the 
global minimum at R = 0, and also due to the friction term, R varies very 
slowly as it approaches the minimum. This explains why in this case the 
model also mimics a cosmological constant. Figure [3] depicts the fraction 
of dimensionless densities f2^ = npi/(3H 2 ) which satisfy the constraint 
fi rad + n ma tt + = 1 where fi ma tt := ^bar + ^dm- The radiation 
contribution, although taken into account, is very small and cannot be 
appreciated from the plots. The current abundances at z — (today) 
match reasonably well the predicted values of the ACDM model and the 
exponential models show an adequate matter domination era. The EOS of 
GDE is plotted in Figure 2] (left panel), and like in other f(R) models [11] , 
it oscillates around the phantom divide value wa = — 1 before reaching 
its asymptotic value as z — > — 1. The total EOS depicted in Figure |4] 
(right panel) shows that at higher z the Universe is dominated by matter 
with wtot ~ 0, and then interpolates to the value u>tot = — 1 in the far 
future. At z = 0, cu to t is similar to the value Wtot ~ 0.75 predicted by 
the ACDM model. Figure [5] (left panel) shows the (modulus) luminous 
distance computed as in and the deceleration parameter (right panel) . 

In these exponential models it is technically difficult to integrate far 
in the past since fun — > exponentially. Since this quantity appears in 
the denominator of Eq. ([5]), it produces large variations that affects the 
precision during the numerical integration. This is something that we had 




Figure 5. (color online) Luminous distance (left panel) compared with 
the Union 2 data pQ. Deceleration parameter (right panel). 



encountered in other f(R) models [TTj . 

The exponential model seems to be consistent with the cosmological 
observations and also with the Solar System [13]. Nevertheless, like other 
f(R) models that look viable as well, a closer examination is required in 
all possible scenarios before considering f(R) theories as a serious threat 
to general relativity. 
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